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Abstract 

Roughly speaking a solitary wave is a solution of a field equation whose energy travels as a 
localised packet and which preserves this localisation in time. A soliton is a solitary wave which 
exhibits some strong form of stability so that it has a particle-like behaviour. In this paper 
we show a new mechanism which might produce solitary waves and solitons for a large class of 
equations, such as the nonlinear Klein-Gordon equation. We show that the existence of these 
kind of solitons, that we have called hylomorphic solitons, depends on a suitable energy/charge 
ratio. We show a variational method that allows to prove the existence of hylomorphic solitons 
and that turns out to be very useful for numerical applications. Moreover we introduce some 
classes of nonlinearities which admit hylomorphic solitons of different shapes and with different 
relations between charge, energy and frequency. 
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A soliton is a solitary wave which exhibits some strong form of stabihty so that it has a 
particle-hke behaviour. 

Today, we know (at least) three mechanism which might produce solitary waves and solitons: 

• Complete integrability; e.g. Korteweg-de Vries equation 

ut + Uxxx + Guux = (KdV) 

• Topological constraints: e.g. Sine-Gordon equation 

utt - + sin u = (SGE) 

• Ratio energy/charge: e.g. the following nonlinear Klein-Gordon equation 

This paper is devoted to the third type of solitons which will be called hylomorphic solitons. 
This class of solitons that are characterised by a suitable energy/charge ratio includes the 
Q-balls, spherically symmetric solutions of INKGi as well as solitary waves which occur in the 
nonhnear Schrodinger equation (see e.g. [10], 0]) and in gauge theories (see e.g. [6], [7], [8]). We 
have chosen the name hylomorphic, which comes from the Greek words ^^hyle^^ = "matter^^ = set 
of particles" and "morphe" = "form" , with the meaning of solitons "giving a suitable form to 
condensed matter" (see Section for definitions and details). 

The aims of this paper are the following: 

• to give the definition of hylomorphic solitons and to set this notion in the literature on 
non-topological solitons; 

• to describe a new variational approach that allows to prove the existence of hylomorphic 
solitons for a large class of nonlinearities for INKGI This variational method turns out to 
be very useful for numerical simulations; 

• to classify the nonlinearities which give hylomorphic solitons: this classification is based 
on the different shapes of the solitons and the different relations between charge, energy 
and frequency. We obtain four classes of nonlinearities and we prove necessary conditions 
for the nonlinear term to belong to a given class. Numerical simulations show the different 
behaviour of these classes quantitatively. 

2 The abstract theory 

2.1 An abstract definition of soliton 

We consider dynamical systems with phase space X described by one or more fields, which 
mathematically are represented by a function 

1- : K" -> y, n>2 

where y is a vector space with norm ||-||y which is called the internal parameters space. We 
assume the system to be deterministic, and denote hy U : R x X X the time evolution map, 
which is assumed to be defined for all t g M. The dynamical system is denoted by {X,U). If 
5*0(2;) G X is the initial condition, the evolution of the system is described by 

5'(i,a;) = Ut'^oix) 
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We assume that X C L2(R",y), namely ^ \\^ {x)\\y dx < +00 for every ^ £ X. Moreover 
for states <E X which satisfy / |a;| da; < +cxd, it is possible to give the notion of 

barycenter of the state as follows 

,it)=Mt) = ^f}^^}^^ (2.1) 

J\\^it,x)\\ydX 

The term solitary wave is usually used for solutions of field equations whose energy is localised 
and the localisation of the energy packet is preserved under the evolution. Using the notion of 
barycenter, we give a formal definition of solitary wave. 

Definition 2.1. A state € X is called solitary wave if for any £ > there exists a radius 
R>0 such that for alH G M 



\\Ut^ix)fydX- f \\Ut^{x)fydX<S 

where Bfl{q^l,{t)) denotes the ball in M" of radius R and centre in qqi{t). 

The solitons are solitary waves characterised by some form of stability. To define them at 
this level of abstractness, we need to recall some well known notions in the theory of dynamical 
systems. 

Definition 2.2. Let X be a metric space and let (X, U) be a dynamical system. An invariant 
set r C AT is called stable, if for any e > there exists (5 > such that if d{^ ,T) < 6 then 
rf(J7t*,r) < s for all t e M. 

Definition 2.3. A state ^'o is called orbitally stable if there exists a finite dimensional manifold 
r C X with e r, such that r is invariant and stable for the dynamical system {X, U). 

The above definition needs some explanation. Since T is invariant, t/t^'o is in T for all 
t. Thus, since F is finite dimensional, the evolution of \l/o is described by a finite number of 
parameters. Thus the dynamical system (F, U) behaves as a point in a finite dimensional phase 
space. By the stability of F, the evolution Ut'^e of a small perturbation ^'^ of might become 
very different from C/(^'o, but Ut^e remains close to F. Thus, the perturbed system appears as 
a finite dimensional system with a small perturbation which depends on an infinite number of 
parameters. 

Definition 2.4. A state ^'o € A',^'o ^ 0, is called soliton if it is a orbitally stable solitary 
wave. 

According to this definition a soliton is a state in which the mass is "concentrated" around 
the barycenter q{t) e M". In general, dim (F) > n and hence, the "state" of a soliton is described 
by n parameters which define its position and other parameters which define its "internal states" . 



2.2 Definition of hylomorphic solitons 

In this section we will expose a general method to prove the existence of non-topological solitons. 
This method leads in a natural way to the definition of hylomorphic solitons. 

We make the following assumptions: (i) there are at least two integrals of motion, the energy 
£ and the hylomorphic charge TL; (ii) the system is invariant for space translations. In some 
models the hylomorphic charge is just the "usual" charge and this fact justifies this name. We 
add the attribute hylomorphic just to recall that H might not be a charge (as in the case of 
the nonlinear Schrodinger equation). Moreover, in many models in quantum field theory, Ti 
represents the expected number of particles, hence for charged particles it is proportional to the 
electric charge. 
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Now, given the set 

Ma = {'^ eX : n{^) = a} 

In order to prove the existence of sohtary waves and sohtons, we follow a method based on the 
following steps: 

• (S-1) prove that the energy £ has a minimum 5*0 on A4^; 

• (S-2) prove that the barycenter (see (|2.ip ) is well defined and that for any e > there 
exists R > such that 

[ \\^{x)\\ldx- [ \\^{x)\\ldx<e Vvfer 

where T is the set of minimisers; 

• (S-3) prove that T is finite dimensional; 

• (S-4) prove that F is stable. 

These steps will be explained in details for the nonlinear Klein-Gordon equation in Section 
The integrals of motion are used in the definition of hylomorphic solitons. Setting 

mo = lim inf -^-I^Xr (2-2) 

where 

^e-{*e^ : <£} (2.3) 

and assuming that mo > 0, we introduce the hylomorphy ratio A {'i') of a state ^ G X, defined 
as 

A (5-) := — ^i^i— (2.4) 

The hylomorphy ratio of a state turns out to be a dimensionless invariant of the motion and an 
important quantity for the characterisation of solitons. 

Definition 2.5. A soliton '^q is called hylomorphic if 

A(*o)<l (2.5) 

In the following we will refer to (j2.5p as the hylomorphy condition. In the study of hylomor- 
phic solitons, the hylomorphy ratio plays a very important role and it is related to a density 
function which we call binding energy density. 

We assume £ and H to be local quantities, namely, given e AT there exist density functions 
Ps,^ {x) and pn,^ {x) G L^(M") such that 



£W = j P£,* (x) dx (2.6) 

ni^) = J PH^^ [x) dx (2.7) 

Then we introduce the binding energy density defined as 

Pb{x) = Pb,^{x) := [mo \pH.^{x)\ - pe,^{x)\^ (2.8) 

The support of pB,^{x) 

E (*) = {x : pB,^ix)^0}. (2.9) 
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is called the condensed matter region since in these points the binding forces prevail, and the 
quantity 

/ P£,^s,ix)dx (2.10) 

"'S(*) 

will be called the condensed matter. The next proposition justifies the choice of the name 
"hylomorphic" , since it shows that hylomorphic solitons "contain condensed matter" . 

Proposition 2.6. //A(^o) < 1 then for all t E M. the support I](C/f5'o) of pB,Ut^{x) is not 
empty. 

Proof. It follows by the simple relations 

J PB,Ut^{x) = J [mo \pH,Ut^ix)\ - P£,Ut^{x)]~^ > 

> mo\n I - £ {Ut^) = mo\n W = 

= 77101^(4-) 1(1- A (^')) >0 

If '^{x) is a finite energy field usually it disperses as time goes on, namely 



□ 



However, if A ('5) < 1 this is not the case. 
Proposition 2.7. IfA{^) < 1 then 



limmf l|t^t^'llL-(R",y) = > 



V) 



Proof. Let A {^) = 1 — a for a given a > 0. We argue indirectly and assume that, for every 
e > 0, there exists t such that 

l|C^f*llL~(R-,y) <e 

namely, C/j^' G where is defined as in (|2.3p . Then, by (|2.2p . if e is sufficiently small 

> ^"0 



and hence 

A(C/f*)>l-|>l-a 
Since A {Uf'^) = A ('J) we get a contradiction. □ 

Thus if A (^P) < 1, by the above propositions the field '5 and the condensed matter will not 
disperse, but will form bumps of matter which eventually might lead to the formation of one or 
more hylomorphic solitons. 
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3 The Nonlinear Klein-Gordon equation 



The nonlinear Klein-Gordon equation (NKG) is given by 

□?A + W'{i>) = (NKG) 
where = - A, ijj{t, x) : M x M" ^ C and : C R satisfies 

W{i;) = Fm) 

for some smooth function F : R+ R. Also we have used the notation 

^From now on, we always assume that 

W(0) = W'{Q) = 

Under very mild assumptions on W, the (NKG) admits spherical solitons. In this section, we 
recall some general facts about (NKG) and study the hylomorphic properties of solitons. In the 
next section we study the problem of existence and stability of solitons, and give a classification 
of the nonlinear terms in (NKG) according to the properties of the solutions it admits. 

We recall the pioneering paper of Rosen [18] and [11], [21], [9]. In Physics, the spherically 
symmetric solitary waves have been called Q-balls by Coleman in [12] and this is the name used 
in the Physics literature. 



3.1 General features of NKG 

In this case, the fields are functions 4* = {ip,ipt) with values in y = and equation (|NKG|) is 
the Euler-Lagrange equation of the action functional 

C{il;uV'4j,'iJj)dxdt (3.1) 



where the Lagrangian density is given by 

C^^l^Ptf-lm'-Wii;). (3.2) 

Moreover the nonlinear Klein-Gordon equation admits a formulation as an infinite dimensional 
Hamiltonian dynamical system. The phase space is given by X = H^{W\C) x L^(M",C) and 
the state of the system at the time t is then defined by a point 5* = ijjt) E X, where ip{t, x) 
is an admissible function for the functional (\3.1\i and t/jt := dtip. Equation (NKG) written as a 
first order system on X takes the form 

a* 



where 




Hence, looking at ipt as the variable conjugate to ^p, the dynamics on X is an infinite dimensional 
Hamiltonian system with the Hamiltonian function given by 

H{^) = ^\^Ptf + ^\Vi;\^ + WW. 
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The nonlinear Klein-Gordon equation is the simplest equation invariant for the Poincare group 
and the action of on X given by 



S^xX3 [9, ^ Tg (*) = (e*V, e^Vt) G X 



(3.3) 



By Noethcr theorem the existence of conservation laws implies the existence of integrals of 
motion. In particular, the time invariance of the Lagrangian implies the preservation of the 
energy E given by 



^l^tl' + i|V7MViy(V') 



dx 



whereas the invariance with respect to the 5^-action (13. 3[) implies the preservation of the hylo- 
morphic charge Ti, given by 



In the analysis of the behaviour of solutions of (NKG), it turns useful to write in polar 
form, namely 

V'(t,x) =w(t,a;)e'^(*'^) (3.4) 

where u{t,x) G K"*" and S{t,x) G M/(27rZ). Letting ut :— dtu, k{t,x) VS{t,x) and uj{t,x) := 
—dtS{t,x), a state S X is uniquely determined by the quadruple {u,ut,iLJ, k). Using these 
variables, the action (|3.ip takes the form 



S{U, Ut,UJ,k)^^ I I Vu|^ + (cj2 - |fc|2) u2 



dx dt 



W{u)dxdt 



and equation (|NKGp becomes 

Uu + (|/.f -u?)u + W'{u) = 

dt {uju^) + V • (ku^) = 
In this form, energy and hylomorphic charge become 



£{u,ut,ui,k) 



\ul + \\Vu\' + \{u;^ + \k\') u^ + W{u) 



dx 



H{u,ut,uj,k) — — / Lou^dx 



(3.5) 
(3.6) 

(3.7) 
(3.8) 



We now describe a possible interpretation of the hylomorphic quantities such as the hylo- 
morphic charge, the hylomorphy ratio etc. The equation (NKG) describes how the density of 
hylomorphic charge moves with the dynamics. Indeed, looking at (NKG) in polar form, one no- 
tices that p.6p is the continuity equation for the density PH,tp of the hylomorphic charge Ti (see 
(|2.7p and p.Sp ). When pn,^ is negative, it can be interpreted as the density of "antiparticles" . 
^From this point of view, equation (|3.6p describes the conservative transport of the hylomorphic 
charge by a particle field and the quantity toq defined in (|2.2p can be considered as the rest 
energy of each particle. 

In this interpretation, TL {^>) represents the total number of particles counted algebraically 
(particles minus antiparticles), \ t}(^) \ represents the average energy of each particle, and the 
hylomorphy ratio A (5*) defined in (|2.4p is a dimensionless quantity which normalise the average 
energy for particle with respect to the rest energy toq- If A {^) > 1, the average energy of each 
particle is larger than the rest energy; if A {^) < 1, the opposite occurs and this fact means that 
particles interact with one another by an attractive force. In this interpretation, the volume 
occupied by the condensed matter defined in p.9p consists of particles tied together. Moreover 
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the quantity pb,<s/ defined in (|2.8p . has the dimension of the energy and represents the binding 
energy of the particles. In fact, consider for example the case of identical particles. When the 
particles are free and at rest, their energy is given by the number of particles times the energy 
of a particle, namely mo|7i (^') |, whereas the energy of the configuration 5" is given by £ (5'). 
Thus the energy necessary to "free" the particles is mo|7i| ~£. Moreover / pB.^{x)dx represents 
the portion of the binding energy which is localised in the support 
Next we will compute the hylomorphy ratio A for (NKG). 

Theorem 3.1. For (NKG) with W of class the hylomorphy ratio defined in I2.4\l takes the 
form 



AW 



Proof. Let us write 



W{s) = -m^s^ + N{s) 



s e 



(3.9) 



where nP = W"{0). We have to prove that the term mp defined in (|2.2p satisfies mo = m. 
Using the polar form for * and ([37l) . ([3^ and (pJj) . we get 





I 




dx 






1/ uju'^dx\ 





> 



> 



J [^w^u^ + im^u^ + N{u)] dx 



J \u}\ v?dx 

Since by the classical Cauchy-Schwarz inequality 



j \lo\ u^dx < (^J ( 



u^dx 



1/2 



Pdx 



1/2 



< 



1 

m 

1 
2m 

1 
2m 



P u^dx 



1/2 



nPu^dx 



1/2 



< 



uPv^dx + j rrpu^dx 
(tj^ + m^) u^dx 



we have that 



^(*) ^ J [k^^u'' + krn''u'^ + N{u)]dx 



J N{u)dx 



f iuP- + vn?) iPdx 



2^ J {up + m^) v?dx 



Then since 7V(u) — o(%P) for u — > 0, it follows that 

mn = lim inf - — ^- — ^— > m 

In order to prove the opposite inequality, consider the family of states "^s.r = {^ur, —imSun), 
where 

1 if |x| < R 

ur{x) = { if |a;| > i?+ 1 

l + if i?< bl <i? + l 
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Then 



< 



e.R) 



dx 



< m + 



m J uj^dx 
1 J\VuR\^dx J N{euR)dx 



< 



2m J uj^dx e^m J uj^dx 



m + 



o{e). 



dx 



where we have again used N{u) = o(m^) for u ^ 0. 



□ 



4 Q-balls 

This section is devoted to the existence of Q-baUs and to the study of their structure. 
4.1 Solitary waves for NKG 

The easiest way to produce sohtary waves of (|NKGP consists in solving the static equation 

-Aw + Ty'(w) = (4.1) 
and making a change of the frame of reference to give a velocity to the wave, setting for example 

tpv{t,x) = 1p^{t,Xi, ...,Xn)^U \ --j==, . . .,Xnj (4.2) 

which is a solution of (jNKGp which represents a bump travelling in the a;i-direction with speed 
V (we consider a normalisation of units of measure so that the speed of light c is equal to 1). 

Unfortunately Derrick, in the very well known paper |13| . has proved that the request W{u) > 
(which is necessary if we want the energy to be a non-negative invariant) implies that equation 
(|4.ip has only the trivial solution. His proof is based on the Derrick-Pohozaev identity, which 
for (|4.ip is given by 

and holds for any finite energy solution u of equation (|4.ip (for details see also [5^ ) . Clearly the 
above equality for n > 3 and W{u) > implies that u = 0. 

However, we can try to prove the existence of solitons for (jNKGp exploiting the possible 
existence of standing waves, namely finite energy solution of (jNKG[) of the form 

Mt,x)^uo{x)e~'^"', uo>0. (4.3) 

Substituting (|4.3p in (jNKGp . one finds that uo{x) is a solution of 

~ Au + W' (u) = Lu^u (4.4) 



and the existence of non-trivial solutions of ()4.4p is not prevented by the Derrick-Pohozaev 
identity, which in this case reads 
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Since the Lagrangian (j3.2p is invariant for the Lorentz group, we can obtain other solutions 
ipv{t, x) just making a Lorentz transformation on it. Namely for n = 3, if we take the velocity 
V = (u,0,0), \v\ < 1, and set 

t' = J {t — VXi) , x[ = J {xi — Vt) , x'2 = X2, = X3 

with 

_ 1 

it turns out that 

■il^^{t,x) := 'ip{t',x') 

is a solution of (jNKGp . 

In particular, given a standing wave ipoit, x) = uo{x)e^'"^°^ , the function il>v{t, x) = V'oC^'j x') 
is a solitary wave which travels with velocity v. Thus, if for example v = (u,0,0) and u(x) — 
u{xi,X2,X3) is any solution of equation (|4.4p . then 

ij^{t,xi,X2,X3)^u{j{xi^vt),X2,X3)e'^^-''-^*\ (4.6) 

is a solution of (jNKGp provided that 

uj — jluq and k = ^ujav 

Notice that (|4.2p is a particular case of (|4.6p when luq — 0. 

Finally we remark that if m is a solution of equation (14. 4p . then by a well known result 
on elliptic equations }14j it has necessarily spherical symmetry. Coleman called the spherically 
symmetric solitary waves of (NKG) Q-balls {12.) and this is the name generally used in the 
Physics literature. 

In order to prove the existence of solitons for (NKG), we are now led to study the existence 
of orbitally stable standing waves. Hence in particular we study existence results of couples 
(uqjUJq) which satisfy equation (|4.4p under some general assumptions on the function W. 

4.2 Existence results for Q-balls 

From now on we make the following assumptions on W: 

• (W-i) (Positivity) W{s) > 0; 

• (W-ii) (Normalisation) 1^(0) = W'{0) = 0; W"{0) = 1; 

• (W-iii) (Hylomorphy) Aq := inf V4 < 1- 

seR+ 2" 

Let us make some remarks on these assumptions. 

(W-i) implies that the energy is positive for any state. Some aspects of the theory remain 
true weakening this assumption. In particular there are results in the case n > 3 also for 

W{s)^ls^~-sP, 2<p<^, (4.7) 
2 p n — 2 

since they are easier to prove with variational methods. However assumption (W-i) is required 
by most of the physical models and it simplifies some theorems. 

(W-ii) is a normalisation condition. In order to have solitary waves it is necessary that 
W"{0) > 0. There are results also in the null- mass case W"{0) — (see e.g. and [2]). 
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However the most interesting situations occur when W"{0) > 0. In the latter case, re-scahng 
the independent variables {t,x), we may assume that W"{0) = 1. 

(W-iii) is the crucial assumption. As we will see this assumption is necessary to have states 
with A (^f) < 1, hence hylomorphic solitons according to Definition 12.51 

Under the previous assumptions, we can write W as 

W{s) = ^s^ + N{s), seM+ (4.8) 

and the hylomorphy condition can be stated by saying that there exists sq G IR+ such that 
-/V(so) < 0. Actually, by our interpretation of (NKG) (see Section [5TT|) . N is the nonlinear term 
which, when negative, produces an "attractive interaction" among the particles. 
In [3] we prove that 

Theorem 4.1 ([3]). // (W-i), (W-ii) and (W-iii) hold then (NKG) m W\ with n>2, admits 
hylomorphic solitons of the form (J^- 



This result is recent in the form given here. In fact only recently it has been proved the 
orbital stability of the standing waves (|4.3p with respect to the standard topology of AT = 
iJi(R",C) X L2(R"^C) and for all the W which satisfy (W-ii), (W-hi) and a condition weaker 
than (W-i) . Nevertheless Theorem 14.11 has a very long history starting with the pioneering 
paper of Rosen [H] . Coleman TT] and Strauss ^21] gave the first rigourous proofs of existence 
of solutions of the type (|4.3p for some particular W, and later necessary and sufficient existence 
conditions have been found by Berestycki and Lions [S]. 

The first orbital stability results for (NKG) are due to Shatah, Grillakis and Strauss [TO], 
[IS] . Namely, they consider the real function 

u i-^ d{uj) := £{u^,u;) + uj\H{u^, cj)| 

where u^^ is the ground state solution of (|4.4p for a fixed to, and £{ui^,lu) and Ti.{u^,Lu) are the 
energy and the hylomorphic charge of the standing wave ipi^{t,x) = u^(x)e^^'^^ . They prove 
that a necessary and sufficient condition for the orbital stability of the standing wave ijj^g is the 
convexity of the map d(uj) in loq. It is difficult to verify theoretically this condition for a general 
W since d{uj) cannot be computed explicitly. In some particular cases the properties of d{uj) 
can be investigated, see for instance [20' where the authors give exact ranges of the frequency uj 
for which they obtain stability and instability of the respective standing waves for the nonlinear 
wave equation with W{ip) as in (|4.7p . Moreover the computation of d{uj) in a neighbourhood of 
a chosen ujq can be extremely difficult due to the fact that the knowledge of is required. 

The proof of Theorem 14.11 follows from some interesting preliminary results, which give the 
idea of the proof. We now sketch the main steps, following the general method exposed in 
Section [^T^ referring to [3J for details. 

The first step is to define the following set 

Y := H\W\R+) X M 
which corresponds to a subset of X through the identification 

Y 3 {u{x),Lu) I — y {u{x),-iuu{x)) eX = H\W\C) x l2(M",C) 

Notice that the standing waves are contained in Y . 

We introduce the functionals energy and charge on the space Y which with abuse of notation 
we still denote by £ and Ti. They take the form 



£{u, ll>) 



^\^u\^ + ^uj^u^ + W{u) 



dx (4.9) 
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H(u,w) = u'^dx (4.10) 

Without loss of generality we restrict ourselves to the case H < and lu > 0. The relation- 
ship between the (NKG) equation and the energy £ on the space Y is given by the following 
proposition 

Proposition 4.2 {'3''). The function tp(t,x) = u(a;)e~"^* is a standing wave for (NKG) if and 
only if{u,u!) is a critical point (with u 0) of the functional £{u, uj) constrained to the manifold 

Ma = {{u,uj)eY : \H{u,uj)\^a}, a^O 

Proof. A point {u,uj) e F is critical for £{u,llj) constrained to the manifold Ma if and only if 
there exists A such that 

au ^^^^^^ 

d£{u^Lj) ^ d'H{u.Lj) 

duj duj 

By definition of energy and charge (equations (|4.9p and (jl.lOp ). equation (14. lip can be written 
as 

-/\u + W'{u) + uj'^u = -2\uju 

Juju'^dx — —Xju^dx. 
Since <t 0, from the second equation A = —uj, and the first becomes equation (|4.4p . □ 

This result gives a simple criterion to obtain standing waves for (NKG). Moreover we are 
interested in standing waves which are orbitally stable (see Definition 12. 3p , and to this aim we 
look for points of minimum of £ on M^- To this aim we use the notion of hylomorphy ratio, 
which on Y takes the form 

. / X ^iu, Lu) I 1 / X . , , „x 

A{u,u; = \ \ = u; + —a{u 4.12 
\n(u,uj)\ 2 2tj 



where 



One immediately gets 



f a\vu\^ + w{u))dx , , 



inf A{u,uj) = Ja{u) (4.14) 



and if Aq is defined as in (W-iii) , then it is proved in 3 that 

Ao = inf a(u) (4.15) 
ueH^ 

To prove the existence of point of minimum for the energy we use 

Proposition 4.3 ([3]). If there exists (u, (D) G Y such that A{u,oj) < 1, there exist points of 
minimum of the energy £{u,uj) constrained to the manifold Ma with a — \Ti.{u,uj)\. 

Equations (|4.14p and (|4.15p imply that assumption (W-iii) is sufficient for the existence of 
{u,uj) with A{u,uj) < 1. Hence step (S-1) of Section [^T^ is completed. 

A point of minimum {uqjLOq) G Y corresponds to the standing wave tpoltjx) = Uo{x)e~^'^°* , 
with uo{x) a spherically symmetric function. Steps (S-2) and (S-3) follow from the fact that for 
isolated points of minimum the set of minimisers F consists of the set 

r := {{uo{x + a)e'\-iujoUo{x + a)e''^)}^^^^ C X 
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which has dimension n + 1. In the foUowing we restrict ourselves to the case of isolated points 
of minimum, which is "generic" in the family of (NKG) equations. 

To finish the proof of Theorem 14. II it remains to prove step (S-4), that is the stability of the 
set r. In 15 we prove that the function 



Vi"^) := -minf 

Ma 



is a Lyapunov function on X. This implies that 

Theorem 4.4 ([3]). //(uo,Wo) is a point of local minimum of the functional £{u, uj) constrained 
to the manifold Ma with a = \Ti.{uQ,LUQ)\, then ipQ{t,x) — Uo{x)e^'"^°* is an orbitally stable 
standing wave. 



Theorem 14. II is a pure existence result and gives no information on the charge and frequency 
of the standing waves. In ,3 we show that Proposition 14.31 implies that, under the assumptions 
(W-i), (W-ii) and (W-iii), there exists a threshold value cto > such that for any cr e ((To,oo) 
there are hylomorphic solitons with hylomorphic charge a. The existence of hylomorphic solitons 
for all hylomorphic charges can be obtained by a stronger version of assumption (W-iii). Let us 
consider the condition 

• (W-iv) (Behaviour at s — 0) N{s) < for s e M+ small enough {N{s) is defined in 

Corollary 4.5 ( 3 ). // (W-i), (W-ii) and (W-iv) hold then for any a 7^ there exists a 
hylomorphic soliton for (NKG) with hylomorphic charge \Ti.\ = a. 

We finish this section by giving a remark which is useful for the numerical approach to the 
existence of hylomorphic solitons for (NKG). We have found solitons as points of minimum for 
the energy £{u,llj) on the manifold with fixed hylomorphic charge. Hence we are studying a 
minimisation problem in two variables with one constraint. It is immediate that this problem 
can be translated into a minimisation problem in one single variable with no constraints. We 
will use as independent variable the functions u e H^. 

If (uo, ujq) is a minimiser of £ (u, co) constrained to with a fixed, then it is also a minimiser 
of A(m,w) constrained to M^- Moreover if {u,oj) G Ma, then 



Hence, letting 



in) 



J dx 



Ka (it) := A (u, LOa (u)) = - 



W{u) 



dx 



we can state Theorem 



in the form 



2 / u'^dx 



(4.16) 



(4.17) 



Theorem 4.6. For any fixed a G if uq{x) is a point of minimum for the functional Ka [u) 
defined in j^.j?] ), then il}o(t,x) = Mo(a;)e~*'^''''"-'* is an orbitally stable standing wave of (NKG) 
with hylomorphic charge \'H\ = a and 0Ja{u) given as in ^773|). 



4.3 Numerical construction of Q-balls 

In this section we provide basic details of our numerical method for constructing hylomorphic 
solitons. The numerical results are shown in the next section, where they illustrate a classifica- 
tion of (NKG) equations originally introduced therein. 

Theorem l4.6l is straightforwardly exploited for the numerical construction of Q-balls: we fix a 
hylomorphic charge a G R"*" and look for points of minimum of the functional Ka{u). Firstly, by 
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I'i'l countour plot 




(a) 



(b) 



Figure 1: (a) Surface plot of |^'| for a 2D Q-ball moving with velocity v = (0.9,0) on the (xi,X2) 
plane (recall that the units of measure are scaled so that the speed of light is c = 1). W given 
by (j4.32p . hylomorphic charge a = 300. (b) Corresponding contour plot (detail), highlighting the 
Lorentz contraction along the direction of motion. 

the classical principle of symmetric criticality (see |16j). we can restrict ourselves to the analysis 
of radially symmetric functions u{x) = u{r (a;)), with r{x) = \x\ in M". We then consider the 
following evolutionary problem: 



in which f denotes a chosen upper bound for the r— domain (discussed below) and lo = ^^(m) as 
in (|4.16p . The evolution of u according to (|4.18p is a gradient flow and therefore a non- increasing 
trend for A(j(t) Ka[u{r,t)) is obtained, well suited as the sought minimisation process. 

The problem (|4.18[) is then discretised by a classical line method; in particular, 2nd order and 
1st order accurate finite differences have been respectively used for space and time discretisation 
(see e.g. [12). The chosen charge a is directly enforced at the n-ih time level (rt = 0, 1, 2, . . . ), 
by evaluating the frequency from the corresponding numerical solution through the discrete 
counterpart of (|4.16p . Moreover, time-advancing is stopped when |ci;"+^ — w^l/w" < e^, 
being a predefined threshold (a relative error on K„ has been considered as well) . The proposed 
method manages to efficiently converge by starting from several initial guesses: not only from 
Gaussian (Q-ball like) profiles but also from discontinuous ones (e.g. u-x {Dr), with Dr contained 
within the chosen r— domain and u suitably set for obtaining the desired <j). Finally, the domain 
extreme f is chosen in such a way that it does not affect the numerical results (an a-posteriori 
check might be necessary: the chosen domain must be large enough to contain the soliton 
support, with some margin). 

We remark that we did not try to implement a shooting method, on purpose. Indeed, 
sign-changing solutions are known to be unstable and, based on the given definition of soliton 
(implying stability), such solutions are of no interest in the present study. Conversely, the 
proposed numerical method guarantees to find a soliton, due to Theorem 14.61 

It is worth remarking that, once defined u(r{x)) by means of the aforementioned strategy, 
it is easy to build moving Q-balls through the transformation (|4.6p : a two-dimensional (2D) 
example is shown in Figure [TJ 



dtu{r, t) ^ - 
u{r, t) = 



a dA^ = Au-W'{u) +Lj'^u in [0,f)xR+ 

on {r — f} X R 



(4.18) 
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4.4 Classification of the nonlinear terms 

In this section we introduce four classes of behaviour for the nonhnear term W, all classes 
admitting hylomorphic solitons. The classification is based on the existence and non-existence 
of solitons with small charge H and with big norm for the (NKG) with in a fixed class. 
Moreover we state some general properties shared by hylomorphic solitons. Without any loss of 
generality we consider the case H < and a; > 0. 

Theorem 4.7 (Admissible frequencies). Let (W-i), (W-ii) and (W-iii) hold and ip(x,t) = 
u(a;)e^"^* be a hylomorphic soliton for (NKG) in R". Set 

Ju^dx 

and 



OJl 



1 1 4r(u) 



where Xq is defined as in (W-iii). Then T{u) < j and there exists 77 > such that the frequency 
w satisfies 

( 1 1 I 4r(u)\ /I , 



Proof. We first show that r<fandw<i + iy^l-^. 

We recall that if ip{t, x) is a hylomorphic soliton for (NKG), then its radial part u{x) satisfies 
equation (|4.4p . Hence using the Dcrrick-Pohozaev identity (|4.5p for u it follows that the energy 
(|4.9|) and the hylomorphy ratio (|4.12p can be written as 

£{u,Lo)^ [ -\\/u\'^ +uj'^u'^dx 
J n 

A/ ^ r(M) 

A(u,Uj) = — ^ + LU. 

nuj 

By imposing that A(u,w) < 1, a simple computation implies the first part of the theorem. 

It remains to prove that uj > lo\. By Proposition 14.151 and (|4.12|) any hylomorphic soliton 
fulfils 

A(U, Lj) = — + - > + - 

Again by imposing that ^ + 7 < 1 we obtain that lu > uji. □ 

Figure [2] shows a "typical" behaviour of the frequency and the hylomorphy ratio of a hylo- 
morphic soliton as the hylomorphic charge varies. The frequency lu and the hylomorphy ratio A 
are decreasing functions of the hylomorphic charge. We remark that as expected we find A < 1 
for all hylomorphic solitons. 

Definition 4.8. Let W satisfy (W-i), (W-ii) and (W-iii). We say that equation (NKG) (or W) 
is of type (a) if there exists ao > such that any hylomorphic soliton ^p{t, x) of the form (|4.3p 
satisfies 

\mt,-)\\L^(R^)>ao VteM 

Theorem 4.9. Assume that W satisfies (W-i), (W-ii) and (W-iii), and write it as in ^■8\ ). If 

N{s) > Vs e (0, ao) 

then equation (NKG) is of type (a). 
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(a) 



(b) 



Figure 2: Frequency (a) and hylomorphy ratio (b) of 2D Q-balls as a function of the hylo- 
morphic charge 7i. W given by (j4.32p : circles associated with the fohowing charge values: 
{5, 10, 20, 30, 40, 50, 60, 70, 80, 90, 100, 200, 300, 400, 500}. 



Proof. Let x) — u{x)e be a hylomorphic soliton for (NKG) and let us assume that 
there exists to such that ||V-'(io: ■)IIl°°(R") ^ '^o- By the assumption on N, this implies that 
N{^P{to,x)) > for aU x e R". 

However, by Proposition 12. 6[ it follows that hylomorphic solitons have non-empty support 
of binding energy density at any time. That is, writing 

1 ,2 1 



P£,^{t,x) = -\Vuf + -[l+uj^]u^ + N{u) 



2 



see equations (|2.8p and (|2.7p . for any t G M there exists a positive measure set of points x E M" 
such that the binding energy density fulfils 

PB,^ ^^v?-\^ |Vu|' - ^ [1 + - ^(w) > (4.19) 

see (I2.10p and Theorem 13.11 which implies mo = 1. Moreover by equation (|4.19p it follows by 
simple computations that for alH £ M there exists a positive measure set of points x G M" such 
that 



< PB,i,{t,x) < 



1 1 2 

ti) — - — -LlI 

2 2 



u^{x) - N{u{x)) < -N{u{x)) 



Letting t — Iq in the previous inequality we get that N{^{to, x)) < for some x, which is a 
contradiction. □ 

Definition 4.10. Let W satisfy (W-i), (W-ii) and (W-iii). We say that equation (NKG) (or 
W) is of type (non-a) if for any ag > there exists a hylomorphic soliton ip(t,x) of the form 
dOl) for which 

\\m-)\\L^(M^)<ao VtGR 

Theorem 4.11. Assume that W satisfies (W-i), (W-ii) and (W-iv), and that Xq defined in 
(W-iii) is positive. Then equation (NKG) is of type (non-a). 

Proof. By Corollarv l4.5l if (W-i), (W-ii) and (W-iv) for W hold, there exist hylomorphic solitons 
of the form (|4.3p for arbitrary small hylomorphic charges. Moreover by Theorem 14.71 if Aq > 
the admissible frequencies of the hylomorphic solitons satisfy 

w > 1 - Vl - Ao > 
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Hence we can find a sequence of hylomorphic solitons ipk = Uk{x)e for (NKG) with liylo- 
morpliic cliarge vanisliing and satisfying 



We recall that the (uk) are radially symmetric functions, hence there exist positive constants C 
and jO, only depending on the dimension n, such that (see [9 ) 

KWI <c ""'",^""^ V|a;|>p (4.20) 

\X\ 2 

Moreover the hylomorphic solitons tpk — Uk{x)e~^'^''* satisfy the Dirichlet problems 

r ~Auk=ioluk-W'{uk) mB{Q,p) = {xeW : \x\ < p} 

{ (4.21) 

[uk{x)=ek for x e 95(0, p) 

where by (|4.20p . {sk) are constants such that Sk ^ as k ^ oo. Now, recalling (|4.8p . without 
loss of generality (see Section 3 in [3]) we can assume that 

\N'{s)\ <cisP + 023"^ with ci,C2>0 and 1< p < g < 2* - 1 (4.22) 

where 2* = -^^2- Hence, since G H^{B{0, p)) and by the assumptions ||Mfe||Hi ^ as fc ^ oo, 
applying the classical Sobolev embedding theorems it follows that the right hand side of (|4.2ip 
satisfies 

f{uk):--{ujl-l)uk-N'{uk)^0 inL' Vse(^l,^^ (4.23) 

where q is given in (j4.22p . Moreover, the classical theory of elliptic regularity (see for example 
[1]) implies that if the functions Uk are solutions of (|4.2ip with f{uk) in for some s, then 

Uk e W^-" and \\uk\\w^.^ < const {Wfkh^ + £fc) (4.24) 

where the constant only depends on the radius p. Now, a classical bootstrap argument applies, 
and it follows that equations (I4.23P and ()4.24p hold for bigger and bigger values of s. We sketch 
the ideas of the first step of the bootstrap. From (|4.23p and (|4.24p it follows that 

/ 2* 

lim ||ufc||vK2.= = VsG 1, — 

k-^oc \ q 

which implies that we can write (j4.23p for s in an interval (1,91) with qi > just applying 
again the classical Sobolev estimates. But then f{uk) belongs to L'^ for all s E (l,gi), then 
(|4.24p hold for all s S {l,qi). This argument can be repeated over and over until equations 
((4:23)) and (14241 hold for all s > 1. 

In particular it follows that (I4.23P and ()4.24p hold for s > and the Sobolev theorems 
imply that 

\Wk\\L^{B{o,p}) < const ||ufc|lvy2,a(B(o,p)) < const (||/fc||L=(i3(o,p)) + Sk) (4.25) 
Hence, putting together (|4.20p and (|4.25p it follows that 

lim ||Mfc||L-(K") = 

and the theorem is proved. □ 
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Definition 4.12. Let W satisfy (W-i), (W-ii) and (W-iii). We say that equation (NKG) (or 
W) is of type (/3) if there exists /3o > such that any hylomorphic sohton ip{t,x) of the form 
(|4.3|) satisfies 

||^(i,-)l!Lo.(R„) </3o VteR 

Theorem 4.13. Assume that W satisfies (W-i), (W-ii) and (W-iii), and write it as in 
If 

N'{s)>0 Vse(/3o,+oo) 
then equation NKG is of type (j3). 

Proof. Let ijj(t,x) — u(a;)e~*'^* be a hylomorphic sohton, then m is a solution of (|4.4p . Setting 
= /^o + I* it is sufficient to prove that v(x) < 0. Let fl := {x : v (x) > 0}. Substituting in 
((441) we get 

\ v^O indn 
Multiplying both sides by v and integrating in fJ, we get 

\S/vf + (1 - cj2)(/3o + v)v + N'{(3o + v)v - w^^l dx> ( {-uj^v^)dx 





since < 1 by Theorem 14. 71 and by assumption on A^'. This implies that u = in il. □ 

Definition 4.14. Let W satisfy (W-i), (W-ii) and (W-iii). We say that equation (NKG) (or 
W) is of type (non-[3) if for any /3o > there exists a hylomorphic soliton iplt^x) of the form 
gSl) for which 

||V'(i,-)llL=o(R„) >/3o VieR 
Theorem 4.15. Assume that W satisfies (W-i), (W-ii) and (W-iii), and that 

W(s) 

lim — ^ = and W{s) > 0, s^O 

s— >oo s 

Then equation (NKG) is of type (non-j3). 

Proof. The argument is by contradiction. Let us assume that there exists a constant C such 
that all hylomorphic solitons have L°° norm bounded by C . Then we can choose a sequence of 
hylomorphic solitons {uk,ujk) with charges Ti.{uk,ujk) diverging and with the properties 

ul dx oo (4.26) 
ll"fe(a;)||L~(K") <C (4.27) 
Let us now consider the sequence of radially symmetric triangle shaped functions Vk defined by 

Vk{r)^\ - -' (4.28) 

1 if |r| > pfc 

with parameters and pk suitable to satisfy / u^dx ~ J vf.dx. For these functions it holds 

We now use ()4.12p and ()4.13p to obtain a contradiction. By the assumptions on W and 
(I4.26P and (I4.27|) . the sequence (ufc) satisfies 

,1 1 fW{uk)dx ^ 1 1 W(C) 
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letting the functions Uk converging to the constant function C. 
On the contrary, by (|4.28p the functions Vk satisfy 



1 1 Ji\Vvk\' + Wivk))dx 1 ^ f WjMk) , 1 
Aiv,,u;,)^-uj, + - j-^^ + + _ 

Hence, letting Mk —>■ oo and pk ^ oo under the condition J uldx = J v'^dx, we get that 

A{vk,ujk) < A{uk,ujk) 
for k sufficiently large. Hence, since T-L(vk,i^k) — 'H{uk,ujk), 

£{vk,ujk) < £{uk,Wk) 

for k sufficiently large, which contradicts the fact that {uk,ujk) are points of absolute minimum 
for the energy. □ 

Given the definitions of "types" , we classify the nonlinear Klein-Gordon equations in terms 
of the following classes 

• (a, 13); if it is of type (a) and type (/?) 

• (a, non-(3); if it is of type (a) and type (non-(3) 

• (non-a, (3)] if it is type (non-a) and of type (/3) 

• (7); if it is of type (non-a) and of type (non-(3) 

In particular the equation (NKG) is of type (7) if for every 70 > there exist hylomorphic 
solitons ij}s{t,x) and ipbit,x) satisfying 

\mt, OIIl-ck") < 70 < WMt, •)IIloo(r") v< e k 

We now give examples of nonlinear terms for all classes, presenting numerical results obtained 
by means of the strategy discussed in Section 14.31 The numerical investigation is carried out in 
a 2D context for simplicity. 

In particular, in Figure [3] we have plotted the radial profiles u{r) of some Q-balls associated 
to the (NKG) equations of the four classes. For each class, the plotted profiles correspond 
to different values of hylomorphic charge. We remark that for (NKG) equations of type (a), 
hylomorphic solitons with arbitrary small hylomorphic charge do not exist. In particular, for 
the choice (|4.29|) . the threshold is approximately at cr = 25. In this case, if the charge is e.g. 
a = 15, the numerical algorithm leads to a radial profile which is not a solution of the Klein- 
Gordon equation in the unbounded domain, but rather a solution of the Dirichlct problem in 
the interval fixed for the numerical analysis. For what concerns hylomorphic solitons with big 
charges, we point out the difference in the shape of the radial profiles for (NKG) equations of 
type {(3) or (non-(3). In the first CclSC, clS expected, there exists a bound for the L°° norm of the 
radial profiles, whereas in the second case the L°° norms increase arbitrarily. Finally in Figure [4] 
we have plotted the energy density (pg,*, solid), the hylomorphic charge density {pn,^, dashed) 
and the corresponding binding energy density {ps,^, dashed-dotted) for some Q-balls of all the 
four classes. 

Example no. 1: W of type (a,/3). The equation 

^tt -Aib+{l + \iP\-a \^\^ + IVI') ^ = 
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is of type (a, (3) provided that a = 2.5. Here 



We have ao = 1 and /3o — 2.5. 

Example no. 2: of type (a, non-(3). The equation 



i^tt - At/) + 



l-IV'l + IVI 



is of type (a, non-13). Here 

M^(V^) = ilog(l-|V| + |^|') 



V3 



arctan 



V3 



(4.29) 



(4.30) 



Example no. 3: W of type (non-a, (3). The equation 

V'tt- AV+ (l-a + ^ = 
is of type (non-a, (3) provided that a E (0,2). Here 



1 ■ '2 a ,_,_,3 1 



(4.31) 



Example no. 4: W of type (7). The equation 



iptt - A^Ij + 



l + l^l 



is of type (7). Here 



w i^p) = 1^1 - log (1 + ivi) = I - 1 + o(ivn 



(4.32) 



We notice that W{s) fulfils the assumption (W-iv). The existence of hylomorphic solitons with 
arbitrary charge is guaranteed by CoroUarv 14.51 
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